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\o • 

^ Recently, a free field construction for correlation functions of the (Z-invariant) eight-vertex model^'^ 

I has been proposed'^ within the algebraic approach to integrable models of statistical mechanics. ^■"'^^ The 

■ free field representation provides explicit formulas for multipoint correlation functions on the infinite lattice. 

00 [ However, the resulting expressions given in terms of a certain series of multiple integrals turn out to be rather 

0^ ■ cumbersome. In this letter we give an explicit expression for the nearest neighbor two-point correlation 

function in terms of a single two-fold integral, and perform some checks. We also discuss the independence 
of the integral representations of the free parameter uq of the vertex-face correspondence entering into the 
Q^l free fields construction.'^ 

1^ ■ Let us briefly recall the notations used in Ref. 3 (see Ref. 1 for a complete definition of the eight-vertex 

model). The fluctuating variables e = ±1 are situated at edges of the square lattice. To each configuration 
of variables £1, 62, £3, £4 ordered around a vertex a local Boltzmann weight R^H^ associated, as it is shown 
\ in Fig. la. The nonzero Boltzmann weights can be parametrized as follows^ 

^ : R{u)++ = R{u)ZZ = a{u) = -ip{u)0,{i^;i^)04i^u;i^)9,{ij^{l-uy,i^), 

R(u)+Z = R{u)Z+ = 6N = -ip(w)04(if;i2f)0i(ifu;i3f)04(if(l-«);i^), 

Riu)^l = Riu)t^ = c{u) = -ip{u)e^{i^;i^)e4ii^u;i^)e,{i^{l-uy,i^), 



Riu)-- = Riu)tt = -d{u) = -ip{u)e, (if ;i2f ) 9^ (if u;i^) 9, (if (1 - u); 



2tr^ 



where 9j (w; r) is the standard jth thcta function with the basic periods 1 and r (Im r > 0) . The normalization 
factor p{u) is irrelevant for correlation functions. 

For definiteness, let us consider the model in the antiferroelectric phase restricting values of the param- 
eters e, r, u to be real numbers in the region e>0, r>l,— l<u<l. For fixed r the parameter e measures 
deviation from the criticality. In the limit e ^ the model has a second order phase transition. In the 'low 
temperature' limit e — > 00 the system falls into one of two ground states (Fig. lb) indexed by i = 0, 1. 

Let Pi'^ be the probability in the thermodynamic limit that the spin in the 'central' edge is fixed to be 
£. The label (i) indicates that spins at the edges situated 'far away' from the origin are the same as in the 
ith ground state so that in the low-temperature limit e cxd the probability is nonvanishing for e — (— )^ 
It has been shown in Ref. 3 that the one-point correlation function 
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(a) (6) 
Fig. 1. (a) definition of the weight matrix; (fe) two degenerate ground states. 
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Fig. 2. Probability Pi*e2(ui,U2): (a) definition; (b) equivalent form obtained by rotating the line vi. 
is recovered from the bosonization procedure. The resulting integral representation 

9i 



1^4(0) J Co 27ri hi{v) 



(2) 



can be reduced to the Baxter-Kelland formula for the spontaneous staggered polarization.^ Here we used 
the notations hj{u) = 9j(u/r;m/er) and 'djiu) = 9j{u;m/e). The integration contour Cq goes over the 
imaginary period of the theta functions (from some complex Vq to Vq + m/e) so that — 1 < Ret; < 0. 

Consider now the inhomogeneous eight-vertex model on the lattice where the spectral parameters in 
two adjacent rows are u\, U2 respectively, while e, r are the same for all sites. ^ Let Pe%2{'U'i,U2) denotes 
the probability of the configuration with two fixed variables as it is shown in Fig. 2a. Here {i) fixes the 
conditions at the infinity so that Pe%2 ~* unless £1 — —£2 = (— )' as e — > cx). The main statement of this 
letter is that the free field construction^ allows one to express the nearest neighbor two-point correlation 
function 

g^'^ (Ul -U2)= ^ £i£2Pii^e2 ^2) 

£l,£2=± 



in terms of a two- fold contour integral as follows: 

92\u)- 



2 — „ ,_. Vl{u) 



MO) 

X "^livi - V2) 



hi{u) jc- 

hi{vi - V2 

hi{vi -V2 + I) 



dvi 
f 1) 



dv2 ^4{Vl +V2 - u) 



hi{Vj) 



(Vj) hA{vj) 



(3) 



The integration contours Ci,2 go over the same imaginary period as Co so that — 1 < Reui < « < Rew2 < 1, 
Ret;2 — Reui < 1. 

Briefly, to get this formula we applied the standard procedure of computing traces of bosonic operators 
over Fock spaces'" and then proceeded as in the one-point function case (see Appendix D of Ref. 3) by 
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applying various identities for theta funetions to provide summation over the SOS variables including infinite 
summation over Fock spaces. We will not go into technical details of this procedure since they are more 
cumbersome than instructive. 

Let us only make a remark on the free parameter uq in the free field representation.^ Although the 
explicit formulas for correlation functions in terms of traces of bosonic fields do contain uo, it is evident from 
the physical grounds that correlation functions must be wo-independent. Since the free field representation is 
based on some assumptions, it is important to prove this statement for the resulting integral representations 
for correlation functions. It can be easily checked that these expressions are nonvanishing double periodic 
meromorphic functions of Uo with two generically incommensurate real periods 2 and 2r, which proves their 
uo-independence. We used the fact that g^^ is independent of uq to Gxuq = U2 + ^ in Eq. (3). In addition, 

to make sure, we also obtained the formula for 92*'' with general uq (which turns out to be more cumbersome) 
and checked numerically that it does not depend on uq. 

To support the validity of the integral representation (3) we compared it with the known results. 
• The partition function differentiation method. The probabilities Piila ^^2) are equal to those shown 
in Fig. 2b because of the Z-invariance. So the correlation function (u) can be calculated as 



(i), s id ,9 d , 9 , , , , 



(4) 



Here k is the partition function per site as a function of the Boltzmann weights,^ 

( ^'im m/2\('\ m/2 —m\( m j_ —m _ _Tn _ „—m\ 

log K{a,b, c, d) = log(c + d)+y ~P >^^~P,^ ^ ,f t ^ ^ 

6 V ' ' ' 6V m(l -p™)(l + a;2™) 

with X = e~'^, p = x"^^, z = e'^^i"^"), and the derivatives are taken at the point characterized by the 
parameters w, e, r according to Eq. (1). 

We checked numerically that the results of (3) and (4) coincide at least up to the fifth decimal digit in 
a wide region of values of e and r. 

• The limiting case r ^ 00. In this limit the Boltzmann weights (1) become those of the six- vertex model 
in the antiferroelectric regime. The integral representation for correlation functions of the antiferroelectric 
six-vertex model is known. ^'^ Analytically the formula (3) gives another integral representation in this limit, 
but numerically the results of integration according to both formulas coincide up to the sixth decimal digit. 

• The limiting case e ^ 0. This is the critical region of the eight-vertex model. Using the Baxter duality 
transformation for the Boltzmann weights^ one can map the model in this region onto the six-vertex model 
in the gapless regime whose correlation functions were found in Ref. 11 (see also Ref. 12). Performing the 
duality transform at the level of correlation functions^" one has to compare our answer with the following 
correlation function 

52^^ (/3i ^ ^ -2 < > (/3i,/32) 

in the notations of Ref. 11 with identification v = 1/r, /3j = yku^. The integral representation for this 
quantity can be written as follows^^ 

„JM,ff, _ „ sh/3 /• dvx f dv2 1 
92 yP) — ~ ^: 



vvjS Jc_ 27ri 7(7+ 



shzv/3 Jq_ 2-k\ Jc^ 27ri sh{vi — (3) sh(w2 — /3) 

S\l{vi-V2) -pr sht/t;j ^ 
sh.v{vi — ■i;2 + Itt) -j--!- sht;,- 

Here the contours C± are from —00 to +00. They are chosen in such a way that P + m (resp. (3) is above 
(resp. below) C+ and /? (resp. /? — 7ri) is above (resp. below) C_. By checking it directly one shows that the 
limit r ^ of Eq. coincides with Eq. (5). 



" We would like to thank Prof. J. H. H. Perk for pointing to us the possibility of this check. 
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• The r = 2 case. Under such specification the eight-vertex model is equivalent to a two of non-interacting 
Ising models.^ In this case —g2^ coincides with the nearest neighbor diagonal correlator of the inhomogeneous 
(Z-invariant) Ising model in the ferromagnetic regime^^ (see also Refs. 2, 14, 15), 



i ^4(0;if)^^(if^;if) 
TT 03(0; if) 01 (i^^^; if) 



where crm,n is the spin variable at the site (m, n) of the square lattice. As we show in the Appendix, Eq. (3) 
reduces to this formula at r = 2. 

Wc hope that a similar integral representation can be obtained starting from the free field representation 
for other multipoint correlation functions, in particular, for two-point functions with separation of 2 lattice 
sites or more. It would be also very interesting to understand whether the elliptic algebra approach proposed 
in Refs. 16 would lead to another bosonization prescription and give a direct procedure of obtaining the 
integral representations of the correlation functions of the eight- vertex model. 
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Appendix 

Let us obtain (6) from (3) in the Ising model case r = 2. The integrand of (3) is antisymmetric with respect 
to the permutation of Vi and V2- This allows one to perform the second integration simply by taking the 
residues at the pole V2 = u. Using the identity 

Mu) hiju) _ MO) h'M hijl) hjju + 1) ^ 
i?i(u) /i4(w) i?'i(0)/i4(0)/i4(l)/ii(w + l) ^^"^^ 

which is valid for r = 2, the resulting expression can be represented in the following form 

_o^i (0)^*1(1) 7^,. ^ jr.A^f dv h4{v + l)h4{v - u + 1) 



gX'{u)=2 ')' , )' J{u), J{u)= — , ^ , — (. (A.l) 

^ ' /i4(0)/i4(l) ^' ^' Jc^2mhx{v + l)hx{v-u + l) ^ ' 

The function J(u) satisfies the following determining properties 

j(« + |)=J(«), J{-u) = J{u), •/(|)=^, 

T/ „N -,, ^ /l4(li)ft.4(0) 

J{u + 2 = -J{u) + ' '^)' 
hi{u)hi{0) 

J{u) is regular on the strip — 2 < Reu < 2, 



which fix it completely to be 



J{u) 



1 02(O,i|)0i(i^</;if) 



2e0i(O,if)0i(i^u;if)' 
Passing to the conjugate module in the coefficient at J{u) in Eq. (A.l) one gets (6). 
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